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Abstract 

We consider the zero-electron-mass limit for the Navier-Stokes-Poisson system in unbounded spa- 
tial domains. Assuming smallness of the viscosity coefficient and ill-prepared initial data, we show 
that the asymptotic limit is represented by the incompressible Navier-Stokes system, with a Brinkman 
damping, in the case when viscosity is proportional to the electron-mass, and by the incompressible 
Euler system provided the viscosity is dominated by the electron mass. The proof is based on the 
RAGE theorem and dispersive estimates for acoustic waves, and on the concept of suitable weak 
solutions for the compressible Navier-Stokes system. 

1 Introduction 

Singular limits arise frequently in the process of model reduction in fluid mechanics. In this paper we 
consider the limit of vanishing ratio electron mass/ion mass in a hydrodynamic model for plasma confined 
to an unbounded spatial domain ft C R 3 . 

1.1 Equations 

For a given (constant) density Ni of positively charged ions, the time evolution of the electron density 
n e = n e (t,x) and the electron velocity u = u(t,x) is governed by the system of equations 

d t n e + div x (n e u) = 0, (1.1) 
m e (d t (n e u) + div x (n e u (g) u)) + V x p(n e ) = div x S(V x u) + n e V x <& - m e ^^, (1.2) 
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A$ = n e - N h (1.3) 

where m e is the ratio of the electron/ions mass, p is the electron pressure, $ is the electric potential, r 
is the relaxation time, and § denotes the viscous stress tensor, 

S(V x u) = n(y x ii + V^u - ^div^ul) , n > 0, (1.4) 

see Anile and Pennisi [3J, Jiingel and Peng [TB], [IS]. Moreover, we suppose the electron velocity satisfies 
the slip boundary conditions 

u • n\ au = 0, [S(V*u)n] x n\ 9Q = 0, (1.5) 
and the boundary is electrically insulated, 

V K $-n| an = 0. (1.6) 
As the underlying spatial domain is unbounded, we also prescribe prescribe the far field behavior: 

u(a;) 0, $(x) ^ as |ac| oo. (1.7) 
Our goal is to study the singular limit and identify the limit problem for m e — > under the condition 

• fi R3 m e , or 

• fi/m e — > as m e — > 0. 

1.2 Ill-prepared initial data 

For m e = e 2 , fi e = /J./s 2 , problem (jl.ll - H~7|) is reminiscent of the low Mach (incompressible) limit of the 
Navier-Stokes system that have been investigated in a number of recent studies, see see the survey papers 
by Danchin [8], Gallagher [16], Masmoudi [22] , and Schochet [25], and the references cited therein. The 
zero-electron-mass limit for the inviscid fluid was treated recetly by Ali et al. [2], see also Chen, Chen 
and Zhang [B]. In the latter case, it is shown that the systems becomes neutral, meaning n e — > Ni, while 
the limit velocity field v satisfies a damped Euler system 

div.v = 0, (1.8) 

9 t v + v-V x v+-v + V x II = 0, (1.9) 
r 

supplemented with the impermeability condition 

v • n| 90 = 0- (1.10) 
In [2], [6j, the authors consider the periodic boundary conditions and well- prepared initial data 
n e (0, ■) = N i + e 2 N , £ , u(0, ■) = u , E , div x u . e = 0. 
In this paper, we focus on the ill-prepared data, specifically, 

n e (0, -)=Ni + eN 0t£ , {N^ £ } e>0 bounded in L 2 n L°°(n), (1.11) 

u(0, •) = u , e , {u 0:£ } e >o bounded in L 2 (ft; R 3 ). (1.12) 

In particular, the gradient part of the velocity field will develop fast oscillations in the asymptotic limit 
e -> 0. 
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1.3 Spatial domain 



In contrast with [2], we consider the physically relevant (unbounded) domains with boundaries. Similarly 
to Farwig, Kozono, and Sohr [12], we focus on the class of uniform C 3 domains of type (a,f3,K). 
Specifically, for each point of x £ d£l, there is a function h £ C 3 (i? 2 ), c 3 (R 2 ) < K, an d 

U a ,p, h = {(y, x 3 ) \h{y) - [3 < x 3 < h(y) + /3, \y\ < a} 

such that, after suitable translation and rotation of the coordinate axes, xq = [0, 0, h(0)] and 

n U aJ3 M = {{y, x 3 ) \h{y) - (3 < x 3 < h(y), \y\ < a}, 

Oil n U at p th = {(y, x 3 ) \x 3 = h(y), \y\ < a}. 

Additional hypotheses imposed on the class of domains are stronger for the inviscid limit so we consider 
the two cases separately. 

1.3.1 Hypotheses in the case of constant viscosity 

Since our method leans essentially on dispersion of acoustic waves, we suppose that 

• the point spectrum of the Neumann Laplacian Ajv in L 2 (f2) is empty, 

in particular, the domain must be unbounded. Although the absence of eigenvalues for the Neumann 
Laplacian represents, in general, a delicate and highly unstable problem (see Davies and Parnovski [10]), 
there are numerous examples of such domains - the whole space i? 3 , the half-space, exterior domains, 
unbounded strips, tube- like domains and waveguides, see D'Ancona and Racke [§]• 

1.3.2 Hypotheses in the case of inviscid limit 

The absence of eigenvalues for the Neumann Laplacian is apparently not sufficient to carry over the 
inviscid limit. We need stronger dispersion provided by the so-called L 1 — L°° estimates well-known for 
the acoustic equation in R 3 , cf. Section \7A\ below. More specifically, we focus on the class of physically 
relevant domains represented by infinite waveguides in the spirit of D'Ancona and Racke [9]. We suppose 
that 

fl£\R 3 , Q = B x R L , L = 1, 2, 3, (1.13) 

where 

B C R 3 ~ L is a smooth bounded domain for L = 1,2, fl = R 3 for L = 3. (1-14) 

Obviously, the domains satisfying (|1.13p . (jl.l4j) belong to the class of uniform C 3 domains of type 
(a, 0, K), and the point spectrum of the Neumann Laplacian is empty. A peculiar feature of the present 
problem is that propagation of acoustic waves is governed by a wave equation of Klein- Gordon type (see 
Section [5]), where dispersion is enhanced by the presence of "damping". In particular, we recover the 
L 1 — L°° estimates even in the case of infinite tubes (L = 1) under the Neumann boundary conditions, 
see Section FTTH 
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1.4 Asymptotic limit 

By analogy with the low Mach number limits, we expect the limit velocity to satisfy the incompressible 
Navier-Stokes system with a Brinkman type damping if fi e = const > 0, and the Euler system (|1.8I - [TTTU|) 
in the inviscid limit fj, £ — > 0. 

In comparison with the low Mach number limit, the main difficulty here is the presence of the extra 
term 

n K - Ni 



1 „ » n c ~ Ni„ A i 
-^n e V x <P = — -V^A" 1 

£ £ 



in the momentum equation (|1.2|) . While the gradient component ji-V x ^> can be easily incorporated into 
the pressure in the limit system, the quantity 



Ni 



n R - N 



£ 

should "disappear" in the course of the limit process e — > 0. To achieve this, the dispersive estimates 
based on the celebrated RAGE theorem will be used. 
Another difficulty lies in the fact that the quantity 



Nt 



n 



n„ - Ni 



is (known to be) only locally integrable; for global analysis, it must be written in the form 



N, 



-V A" 1 



Ni 



(div x (V x $ ® V x 4) - iv^V^I 2 



meaning as an element of the dual space W -11 . 

Last but not least, we point out that the analysis of the inviscid limit leans heavily on the fact 
that the propagation of acoustic waves is governed by the Klein-Gordon wave equation, yielding effective 
dispersion on the waveguide like domains specified in Section 11.3.21 

The paper is organized as follows. In Section [2l we introduce the concept of suitable weak solution 
to system (jl.ll - [T7f|) that proved to be very convenient for studying the inviscid limits, cf. [13] . Section 
[3] contains the main results. In Sectional we summarize the uniform bounds independent of the scaling 
parameter e. Section[5]is devoted to the acoustic equation and the resulting dispersive estimates. Finally, 
in Section we show convergence toward the incompressible Navier-Stokes system in the case of non- 
degenerate viscosity, while Section [7] completes the proof of the inviscid limit. 

2 Suitable weak solutions 



Motivated by the genera theory developed by Ruggeri and Trovato [23], we assume that the electron 
pressure p satisfies 



p G C[0, oo) n C 2 (0, oo), p(0) = 0, p'(n) > for all n > 0, lim 



p'{n) 



Next, we introduce the standard Helmholtz decomposition of a vector field v, 

v = H[vl+H ± [vl, 



Poo > 0. 



(2.1) 
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with 

H 1 " = V x *, A* = div^v in O, (V^ - v) • n\ dQ = 0. 
As shown by Farwig, Kozono and Sohr |12) . the linear operator 

H is bounded in L 2 n L q (fl; R 3 ) for 2 < g < oo, and in L 2 + L q (fl) forl<g<2 

as soon as fl is a C 2 — domain of type (a,f3,K) introduced in Section [1731 Moreover, the norm of H in 
the aforementioned spaces depends solely on the parameters (a,f3,K). As a matter of fact, the domains 
considered in the present paper belong to the higher regularity class C 3 for several technical reasons that 
will become clear in the course of the proof of the main results. 

Following |15) we say that a trio n e , u, $ is a suitable weak solution to system (|l.ll - [l~7l) in (0, T) x CI, 
supplement with the initial conditions n e (0, •) = no, u(0, •) = Uo if: 

• the functions n e , u, $ belong to the regularity class 

n e >0, n e — N t e L°°(0, T; L 5 ^ 3 + L 2 (fi)), 

J) ( Jle )ei 1 (0,r;Lj JC (O)) ) 

u G L 2 (0, T; W ia {Vl- R 3 )), u • n\ dn = 0, 
V. T $ G L°°(0,r;£ 2 (O;i? 3 )); 

• equation of continuity (jl.ip is satisfied in the sense of renormalized solutions (see DiPerna and 
P.-L.Lions p]), 



T 

o Jn 



((b(n e ) + n e ) dtp + (b(n e ) +n e )u- V x tp + (b(n e ) - b'(n e )n e ^j div^u^ da; dt (2.2) 



(6(no) + n )) <p(0, •) da; 
for any test function tp £ C c °°([0, T) x H) and any b G C°°[0, oo), 6' G C c °°[0, oo); 



• momentum equation (jl.2p . together with the slip boundary condition (|1.5p . is satisfied in a weak 
sense, 

J J (jn e n e u ■ dt<p + m e n e (u (g> u) : \7 x tp + p(n e )div x <p^ dx dt (2-3) 

= / / (S(Va;u) : V x tp - n e V x $ ■ tp + H^£ u . ^ j da; dt— n u ■ tp(0, •) da; 
Jo Jn ^ T ' Jn 

for any test function e C£°([0,T) x fl;R 3 ), tp ■ n|[ 0iT )xao = 0; 

• the electric potential $ is given by formula 



V*$(s, •) = V x $a - I H x [n e u] dt, 
Jo 



(2.4) 



where 

A$ = n — Ni in f}, V x $ ■ n| an = 0; (2.5) 
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• the relative entropy inequality 

J Qm e n e |u - U| 2 + E(n e , r) + ^V^ 2 ) (s, •) dx (2.6) 

+ / / [S(V B u) -S(V K U)] : Vx(u-U) da; dt + f f — n e \u - U| 2 dx dt 
Jo Jn Jo Jn T 

< J^(^m e n \u ~V(0,-)\ 2 + E(n ,r{Q,-)) + ^\\7 x ^ \ 2 ^j dx + jf ft(n e ,u,r,U) di 

holds for a. a. s £ [0, T] and all test functions r, U such that 

r - Ni G C c °°([0,T] x H), r > 0, U € C c °°([0,T] x H; i? 3 ), U • n\ dn = 0, 

where 

K (n e , u, r, U) = J m e n e (d t V + uV K U) ■ (U - u) da; (2.7) 
+ / §(V K U) : V X (U - u) dai + / — n e U • (U - u) da; - / n e V K $ • U da; 
+ J ((r - n e )d t P(r) + V x P(r) • (rU - n e u) - div x u(n e (p(n e ) - P(r)) - £(n B ,r))) dx, 



with 

E{n e ,r) = H(n e ) - H'(r)(n e - r) - H{r). 



and 



P = If, H(n) = n /" ^jUs. 

ii s 



It can be deduced from (|2.4I) that 

Va,$(i, -)Vx¥> da; = / (iVi - n e )(t, da; for any t G (0, T) and all test functions ip G C™(Ti), (2. 
Jn 



f 



meaning $ is a (strong) solution of Poisson equation (|1.3j) . In particular, by virtue of the standard (local) 
elliptic regularity, 

•) G W 2 ' 5/3 (K) for any compact K C ft. (2.9) 

The existence of global-in-time suitable weak solutions to the compressible Navier-Stokes system in a 
bounded spatial domain and the no-slip boundary conditions was proved in [151 Theorem 3.1] with the 
help of the approximation scheme introduced in [13] . Adaptation of the method to the present problem 
requires only straightforward modifications. The main advantage of working directly with suitable weak 
solutions is that the relative entropy inequality (12.6[) already implicitly includes the stability estimates 
necessary to perform the inviscid limit. 
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3 Main results 



We start by introducing the scaled system. To simplify notation, we set m e = e 2 , n e = n £ , u = u £ 
$ = $ £ , and Ni — n - a positive constant. Accordingly, system of equations (jl.ll - [TT3l) reads 



dtn £ + div x (n e u e ) = 

d t (n £ u £ ) + div x (n £ u £ ® u £ ) + \\7 x p(n £ ) = div x S £ (S7 x u £ ) + \n £ \7 x $ £ - -n £ u £ 

A$ e = n £ — n, 

with the viscous stress 

S e (Va,Ue) = Me (V x U e + V^.U e - ^div^Ujl) , \l £ > 0. 

System (|3.1I - !3-3[) is supplemented with the boundary conditions 

u e ■ n\ dn = 0, [S e (V x u e )n] x n|ao = 0, 
V,$ £ • n|an = 0, 

and 

u £ (x) -> 0, $ £ (x) -> as |x| -> oo. 

3.1 Ill-prepared initial data 

Taking 

r = n, U = 

as test functions in the relative entropy inequality (|2.6p we obtain 



7: n e\u £ \ 
Q \ 1 



H(n £ ) - H'{n){n £ - n) - H(n) 



^|V«3 8 | a )(.s.-)d.r 



o ./n 



Va;U e + Vj,u e - -div K u £ 



|u F r da; dt 



where 
and 



n £ = n 0: e, u e (0, •) = u , e 
$o, £ = A^ x [n , e -n] 



(3.1) 

(3.2) 
(3.3) 

(3.4) 

(3.5) 
(3.6) 

(3.7) 



(3.8) 



(3.9) 
(3.10) 



are the initial data. 

Consequently, the initial data must be chosen in such a way that the expression on the right-hand 
side of (|3.8p remains bounded for e — > 0. Accordingly, we suppose that 

%(0,') = %=S + £% {N , £ } £> o bounded in L 2 n (3.11) 

u £ (0, ■) = u .e, {u 0ie } £>0 bounded in L 2 (f>; i? 3 ). (3.12) 

Moreover, the functions JVo e must be taken so that $o,s = £ ^-n 1 Wo,e] satisfy 



{V x $o, E } e >o bounded in L 2 (Q;R 3 



(3.13) 
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3.2 Asymptotic limit for positive viscosity coefficients 

Our first result concerns the asymptotic limit in the case \x e = \x > 0. 

Theorem 3.1 Let Cl c R 3 be an (unbounded) C 3 — domain of type (a,(3,K) specified in Section ] 1.3\ and 
such that the point spectrum of the Neumann Laplacian An in L 2 (Cl) is empty. Suppose that the viscosity 
coefficient fi E = /j > is independent of e and that the pressure p satisfies t2.1]) . Let {n e , u e , $ £ } £ >o 
be a sequence of suitable weak solutions to the scaled system I13.1\ - \3.T\j , emanating from the initial data 
satisfying iTOTl - EH^) . 
Then 

ess sup \\n £ -n\\ L 2 +L 5/3 {n) < ec, (3.14) 
te(o,T) 

and, at least for a suitable subsequence, 

Uo.e — > Uq weakly in L 2 (Cl; R 3 ), 

u e -> U weakly in L 2 (0, T; W 1,2 (Cl; R 3 )) and (strongly) in L 2 ((0, T) x K; R 3 ) (3.15) 

for any compact K G Cl, where U is a weak solution to the incompressible (damped) Navier-Stokes system 
in (0,T) x CI, 

div x U = 0, (3.16) 
n[d t \5 + U ■ V X U) + V x n = div x S(V x U) - -U, (3.17) 



with 
and 



U • n\ 9n = 0, [S(V a U)n] X n\ dn = 0, (3.18) 
U(0,-) = H[U o ]. (3.19) 



Remark [3ll Momentum equation \3.17\ ), together with the slip boundary conditions 13.18\) and the 
initial condition \3.19\) . are understood in the weak sense, specifically, the integral identity 



/o Jo. N ' "'o Jn 

for any test function (p € C£°([0, T) x Cl; R 3 ), div x tp = 0. tp ■ n\gn = 0. 



/ [ n(\J-d t (p + (U®U) : V x f) dx dt = / j (s(V x U) : V x <p+ -U- ip) dx dt- j U • <p(0, •) da; 

(3.20) 



3.3 Inviscid limit 

Our second result concerns the case of vanishing viscosity coefficient (j, e \ 0. In this case, the limit 
velocity field is expected to satisfy the incompressible Euler system (|1.8I - 11.101) . As is well-known, this 
system possesses a local-in-time solution 

veCflO.U;^^!])), V X U e C([0,T m , x );W k -^ 2 (Cl)) (3.21) 

provided 

v(0, = vo G Vl^< 2 (r!), k > 5/2, div.vo = 0, v • n| 9n = 0, 
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and provided ft = R 3 , f2 is a half-space, or f2 is an (exterior) domain with compact boundary. The life- 
span T max depends solely on ||vo||vK fc . 2 (n;_R, 3 )i see Alazard pQ, Isozaki [17], Secchi [26], among others. As a 
matter of fact, the damping term -v in (|1.9j) may extend the life-span of regular solutions, in particular 
if the initial data are small in comparison with 1/r. In a very interesting recent paper, Chae [5] showed 
that a smooth solution of (|1.8l - ll.l0|) exists globally in time provided — R 3 , and r < T^ax) where T®^ 
is the life span of the regular solution of the undamped Euler system emanating from the same initial 
data. 



Theorem 3.2 In addition to hypotheses of Theorem \3.1\ assume that Q C R 3 is an infinite waveguide 
specified in Section \l.3.2\ Moreover, we suppose that 

He \ 0, 

and that the initial data satisfy 

7Vo, e — > Nq (strongly) in L 2 (f2), Uo. e — > Uo (strongly) in L 2 (f2;i? 3 )) 

as e — > 0, where 

v = H[u ]e W k ' 2 (n;R 3 ), fc>5/2. 

Moreover, suppose that the damped Euler system i 1 . 8\ - \ 1 . RJ\) . with the initial datum vo, possesses a regular 
solution v defined on a time interval [0,T max ) satisfying \3.21\ ). 
Then 

«llz=+£5/3(n) < ec, (3.22) 



ess sup \\n e 
te(o.T) 



and 



ess sup 

*6(0,Ti oc 



H 



-tu 



L 2 (0;_R 3 ) 



^u e -> v in L q ({), T\ oc ; L 2 (K; R 3 )) for any compact K C O, 1 < q < oo, 



/or <m?/ T loc < T max , Ti oc < T. 



Remark [3]2 The proof of Theorem \3.2\ leans essentially on the L 1 — L°° bounds for acoustic waves 
established in Section \7.1\ Thus the conclusion of Theorem \3.2\ remains valid as soon as these bounds are 
available. Note that Isozaki Ji7| / established similar estimates on exterior domains in R 3 . 

The rest of the paper is devoted to the proof of Theorems 13. H 13.21 



4 Uniform bounds 

For the ill-prepared initial data, all desired uniform bounds follow from the energy inequality (|3.8[) . 
Introducing the essential and residual parts of a quantity h, 

h [^]css ~l~ [^]res 
[ft] bbb = X( n e)h, [ft] res = ft - [ft ess ], 



!) 



where 

X G Cc°(0, oo), 1 < x < 1j x( n ) = 1 f« i! belonging to an open neighborhood of n in (0, oo), 

we get the following list of estimates 

ess sup \\^/n^n s (t,-)\\ L 2^ n . R 3) < c, 
te(o,r) 



ess sup \\V x ^ E (t,-)\\ L 2 {n , R3) < ec : 
te(o,T) 



ess sup 

te(o,T) 



(V) 



< c, ess sup / [n e ]rlg(t, •) dx < e 2 c, 
L2(n) te(Q,T)Jn 



ess sup / [l] rcs (i, ■) da; < e 2 c, 
te(o : T)Jo 



and 



, 2 

V x u E + Vo;U e - -div x u e 



da; dt < 



(4.1) 
(4.2) 

(4.3) 

(4.4) 

(4.5) 



>0 JQ 

where all generic constants are independent of e. 

Estimates (|4.3I - 14.5[) can be combined to deduce a bound on the velocity field in the Sobolev space 
L 2 (0,T; W 1,2 (Vl] R 3 )) that is relevant in the proof of Theorem 13.11 To this end, we report the following 
version of Korn's inequality that may be of independent interest. 

Proposition 4.1 Let fl C R 3 be a C 2 -uniform domain of type (a, (3, K) introduced in Section ] 1.31 
Then there exists 5 > 0, depending solely on the parameters (a,/3,K), such that 



V x w + \7 x w — -div^wl 



Il w llwi.=(n;fl3) < c{a,P,8,K) 
for any measurable set V, \V\ < 5, and for all w £ W 1 ' 2 ^; R 3 ) 



2 

4 


/ |w| 2 dx ) 


L 2 (ri;i? 3 x 3 ) 


Jn\v J 



(4.6) 



Proof: 

In view of the standard decomposition technique and partition of unity, it is enough to show (|4.6[) on 
each set 

U a,p,K = {(y> x 3) \Kv) - P < x 3 < h(y), \y\ < a}. 
Revoking the result [H Proposition 4.1], we have 



H w lliU 2 (£/-;fl3) ^ c(a,/3,K) 



V x w + V x w div^wl 



L 2 (U-;R 3x3 ) JU 



I 2 dx 



(4.7) 



As a matter of fact, the constant c in (|4.7j) depends only the Lipschitz constant of the function h and 
width of U~ given in terms of a, f3. 
Furthermore, we have 

\V-p, K \>2S>0 
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for a certain S(a, /3, K) > 0. In particular, 

| \V\ > S for any measurable set V, \V\ < S. 
Now, arguing by contradiction, we construct sequences 

{hn}%Li, HMc 2 (fi 2 ) < K, U~ = x 3 ) h n (y) - f3 <x 3 < h n (y), \y\ < a|, 

h n ^hm C\{\y\ < a), 

\V n \<6, 



and 

such that 



{w n }^L l5 llWnH^i. 2( L r-. fl .3) 



Va;W„ + V^W„ - -div^Wnl 



|w„ I 2 da; — > as n — > oo. 



£ 2 (t/-;.R3x 3 ) ./E/-\V„ / 

Because the domains U~ are uniformly Lipschitz, we can extend w„ as w„ on the cylinder 

U = {(y, x 3 ) | |ar 3 | < 2 + (3 + aK, \y\ < a} 

in such a way that 

l|w„|| W i,2 (C ;. H 3) < c(a,f3,K)\\vf n \\ wl , 2(u -. R3) . 
Since W 1,2 (U ; i? 3 ) is compactly embedded into L 2 (U ; R 3 ), we may use (|4.7[) to deduce that 

w„ -> w in W 1,2 (?7;i2 3 ), 

where 

V x w + V*w- ^div a wl = 0, w^Ointheset |(y,x 3 ) h(y) - P < x 3 < h(y) , |y| < aj. (4 
On the other hand, 

/ V\v; l w »| 2 dx ^ f X|w| 2 dx = 
Ju Ju 



where 



u~\v„ 



X weakly-(*) in x > 0, / x dx > 0. 



(4.9) 



However, relation (|4.8|) implies that w is a (nonzero) conformal Killing vector (see Reshetnyak [23] ) 
vanishing, by virtue (I4.9p . on a set of positive measure, which is impossible. 

Q.E.D. 

Thus, finally, taking V the "residual set", V = supp[l] rcs we may combine the estimates (|3.8j) . (|4.4j) . 
and (|4.5j) with Proposition 14. II to conclude that 



(4.10) 



Me / IKIIwi,2(fW) d * ^ 
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5 Acoustic equation 

As already pointed out, the essential piece of information necessary to carry out the asymptotic limit is 
contained in the oscillatory component of the velocity field responsible for propagation of acoustic waves. 
Introducing new variables 

n £ — n ^ r 
N £ = , V E = n £ u £ 

e 



we can formally rewrite system (|3.ip . (|3.2[) in the form 

Ed t N s + div x V e = 0, (5.1) 

ed t V e +p'(n)V x N e + -V x $ e (5.2) 

= £div x S(V x u e ) - ediv x (n e u e ® u e ) - ^n £ u £ - iV e V x $ e - -V ' x (p(n e ) -p'(n)(n £ -n) -p(nfj, 
supplemented with the boundary condition 

V £ • n| an = 0. 

System (|5.1j) . (|5.2[) is usually called acoustic equation, see Lighthill [21]. Its (rigorous) weak formula- 
tion reads 

/ / (eN e dt<p + V E • V x <p) dx dt = -e [ N Q , e tp{0, da; (5.3) 
for any 93 € C£°([0,T) x O), and 

y" y (eV e • <9 t V^ + p'(n)N e A<p - riN e <p) dx dt (5.4) 
= e [ T I S £ (V,u £ ) : dxdt-eT ( ( n £ u £ u £ + ^) ~ ^ (n)(n £ - n) - P (n) A . y ^ ^ ^ 



Jn Jo Jn 

T 

+e 



y y (^~ n £ u £ ' V ^ + ^e v x ' V ^ dx dt-e J n , 6 u , e • W x (p(0, •) da; 

for any ip £ C£°(f2), V^y • n|an = 0. Moreover, we rewrite 

y f ^ • v x ip dx = i y a^v*^ . d^ 

= ^y f~V x |V x $ e | 2 • V x <p + div x (V x $ e <8> V x $ e ) • V^j dx 

= y Q|V B $ e | a • + (V x $ e <g V x $ £ ) : dx. 

Furthermore, it follows directly from the uniform bounds established in (|3.8I - 14.51) that (|5.4|l can be 
written as 



y y (eV e -d t V x <p + p'(n)N e Ap--fiN e (p) dx dt (5.5) 
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: J J (Gl : V 2 x tp + G 2 : V%tp + ^n e u E • V^) dx dt - e J n , e u , e • V x tp(0, ■) dx 
for any tp G C£°(fi), V x tp ■ n\on — 0, where 

{Gl} e >o is bounded in L°°(0, T; L^fi; i? 3x3 )), (5.6) 

and 

{<G 2 } e>0 is bounded in L 2 (0, T; L 2 (f2; i? 3x3 )). (5.7) 

5.1 Neumann Laplacian 

At this stage, it is convenient to rewrite the acoustic system (|5 . 3[) . (|5.5[) in terms of a single self-adjoint 
operator A in £ 2 (D), specifically, 

A = -p'(n)A N +nl, 

with 

V(A) = {we W ia {Vt) 

I (p'(n)V x w ■ V x tp + nwtp) dx = / gtp dx for any tp £ C^(Q) for a certain g e L 2 (f2) >. 
Given the regularity of the boundary <9f2, it can be shown that 

V{A) = | to G T^ 2 ' 2 (fi) V.w • n| afi = j, -4[w] = -p'(S)Aw + nw. (5.8) 
Furthermore, since fi is of uniform C 3 -class, the classical elliptic theory yields 

v(A 2 ) c G 2+v n w 2 '°°(Ti) (5.9) 

for a certain v > 0. We remark that all we need is only uniform C 2+l/ — regularity of the boundary instead 
ofC 3 . 

5.2 Acoustic equation - abstract formulation 



£ 



In view of (|5.8[) . (|5.9[) . and the uniform bounds established in ()5.6|) . (|5.7p . the acoustic equation ()5.3|) . 
(15.51) can be written in a concise form: 



/ / (eN £ d^ + V £ • V x <p) dx dt = -e f AT o , E ^(0, •) da; (5.10) 

JO JQ V y «/n 

for any tp G C£°([0,T) x O), and 

/ / (eV e • d t V xV ? - iV e ^M) dxdt = e [ [ F e A 2 [p] dx dt - s f n , e u Q . s • V x¥ >(0, •) da; (5.11) 
Jo Jn ^ ' Jo Jn Jn 

for any tp G C x ([0, T); X>(-A 2 )), with 

{F e } s>0 bounded in L 2 (0, T; L 2 (n)). (5.12) 
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Thus, using the standard variation-of-constants formula, we obtain 



N K 



- ( exp (iVA- 



exp WA 



t - s 



exp ( — iVA- 



Nq.e Fr-^O.e 

VA 



(5.13) 



exp — iv A 



t - s 



ds, 



Z E = - exp iVA- Z , e - ivC4[JV 0lE ] + exp -iy/A- Z , e + ivC4[iV , e ] 



-t 



exp iv A 



t- s 



exp -WA 



e 

t - s 



(5.14) 



[A 2 [F e ]} ds, 



(5.15) 



(5.16) 



where Z s is interpreted as 

/ Z e G(A)[<p] dx = - I V e • V*(G(.A)k/>]) da; for any G € C~(ra,oo), p € C c °°(fi), 
in particular, 

/ Z 0lB G?(^)M dx = - [ n , £ u ,e • V x {G(A)[<p\) dx for any G € C c °°(n, oo), p e C c °°(n). 
Jn Jn 

Note that the spectrum of the operator A is the half-line [n, oo). 
5.3 Application of RAGE theorem 

With the explicit formulas f|5 . 1 3[) . (|5.14l) at hand, we are ready to show local energy decay for N £ and 
the acoustic waves represented by the gradient component H x [V e ]. To this end, we employ the following 
version of the celebrated RAGE theorem, see Cycon et al. |jj Theorem 5.8]: 

Theorem 5.1 Let H be a Hilbert space, A : T>(A) C H — >• H a self-adjoint operator, C : H —> H a 
compact operator, and P c the orthogonal projection onto the space of continuity H c of A, specifically, 

H = H c cl#|span{u> £ H \ w an eigenvector of 
- f exp (-it A) CP C exp (1L4) di 

T JO 



Then 



as t — } oo. 



(5.17) 



C(H) 



We apply Theorem O to H = L 2 {Q), A = -</A, C = x 2 G{A), with X e C c °°(rj), x > 0. In 
accordance with hypotheses of Theorem 13.11 the point spectrum of A is empty, and we deduce that 



X G(A)exp(iVA-)[X] 



L 2 (n) 



dt<w( e )||X||i a 



(n)> 



(5.18) 



where lj(e) — > as e — > 0. In particular, going back to (|5.13[) . f|5 . 14[) we may infer that 

tt-> I N e (t,-)G(A)[<p] dx] -> in L 2 (0,T), 
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and, similarly, 

{t^J V e >V m {G(A)[<p])dx\-Kim. L 2 (0,T) 
as e — > oo for any G G C£°(n, oo), </? € C£°(f2). Thus, by means of density argument, 

N e -> in L 9 (0,T; (L 2 + L 5/3 ) wcak - W (r!)) for any 1 < <? < oo, (5.19) 

while 

jt h->- y" n £ u £ • H x [<^] dxj ^ in L q {0,T) for any 1 < q < oo, ^ <= C™(£1). (5.20) 

6 Compactness of the solenoidal part - proof of Theorem 13.11 

In this section, we complete the proof of Theorem 13.11 To begin, we remark that relation (|3.14[) follows 
directly from (|4.3p . while (|4.10p implies that 

u e -> U weakly in L 2 (0, T; W h2 (Q; R 3 ), (6.1) 

at least for a subsequence as the case may be. Moreover, the vector field U is solenoidal and satisfies the 
impermeability condition U • ii\qq = 0. 

Next, the uniform bounds (|4.2[) . (I4.3[) , together with the standard elliptic theory, yield 



ess sup 

te(o,T) 



V, 



< c(if) for any compact if C fi, (6.2) 



which, combined with (|5.19[) . yields 

J J ^7" ^ (~e~) ' V dX dt ° aS £ ° f ° r ^ € ^([O^) x »;i? 3 ). (6.3) 

Taking = C£°((0,T) x O; i? 3 ), div^ = 0, as a test function in the momentum equation (|2.3p and 
making use of (|6.3[) . we deduce that 



{H[n £ u £ ]} £>0 is precompact in C weak _ w ([0,T]; L 2 + L 5 / 4 (f>; R 3 )). (6.4) 
Indeed the only quantity term reads 

1 „ _ ?! » 7T„ $e 

e l e z e e 



where the former term is a gradient, while the latter satisfies 

Putting together (|5.20[) . (|6.1I) . (|6.4I) . with (|3 . 14[) . we deduce the desired conclusion 

u £ -> U (strongly) in L 2 ((0, T) x K; R 3 ) for any compact K C O. (6.5) 



With relations (|6.3p . (|6 . 5|) at hand, it is not difficult to perform the limit e — > in the weak formulation 
of momentum equation (|3 . 2|) to obtain (|3.20|) . 
We have proved Theorem 13. II 
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7 Zero viscosity limit - proof of Theorem 13.21 



Our ultimate goal is to prove Theorem 13.21 The basic tool here is the relative entropy inequality (|2 .6[) 
satisfied by the suitable weak solutions. Taking n e = n £ , u = u £ , r = n in the rescaled variant of (|2 .6[) 
we obtain 

J a f^K - U| 2 + ^E(n e ,n) + ±V X (^j ) (s, •) dx (7.1) 

+ / S / [§e(V x u £ ) - S £ (V X U)] : V K (u £ - U) dx dt + f f -n £ \u £ - U| 2 dx dt 
Jo Jn Jo Jn T 

- St (^ na ^ u °> £ ~ v °> £ \ 2+ ^ E ( n ^ £ ^ + l Vx (% £ ) ) dx+ S 7? -( n£ ' UE ' 7T ' u ) a< ' 

n (n e , u £ , n, U) = J n e (d t V + u £ V x u) • (U - u £ ) (7.2) 

f § £ (V X U) : V K (U - u e ) dx + f -n e V • (U - u e ) dx - \ [ n £ V x <P £ • U dx 
Jo. Jn T £ Jn 

l(n e (p(n e ) - P(n)) - E(n e ,nj) dx. 



with 



4 / <liv,U( 



e Jn 



Furthermore, we take 



U = U £ ,a = v + V^.a, 



where v is the (unique) solution of the damped Euler system (jl. 8111. 10")) , emanating from the initial data 
vq = H[uq], and V x ^e,s mimicks the oscillatory part of the velocity field. Specifically, we take 



ed t s e .8 + A* £j<5 = 0, 
ed t V x ^ £ j + p'(n)V x s St5 - nV x A~ 1 s £tS + -V.f^ = 0, V,*^ • n| a n = 0, 



(7.3) 
(7.4) 



which is nothing other than a slightly modified homogeneous part of the acoustic system (15.101) . (|5.11[) . 
The initial data are taken in the form 



S £ ,S(0, ■) = i[No,e}s, *e,s(°> = [*0,J 3 , with V.fo, = H X [u , e ], 

n 

where the brackets [-]g denote a suitable regularization operator specified in Section I7TT1 below. 
Keeping (|7.3I - 17.41) in mind, we can rewrite the remainder (|7.2I) in the form 



(7.5) 



n(n e ,u s ,n,U EtS ) (7.6) 

/ n £ (d t v + u £ V x v + -v) • (U e ,5 - u e ) dx + [ (n £ d t V x ^ £ .s • v + n £ u £ ■ V 2 x ^ £ .s ■ (U £> 5 - u e )) da; 
isi v T ' Jn 



SeCV^Ue.5) : V x (U e>5 - U e ) dx - / A* £ ,5 

n Jn 



P{n £ ) — P(n) n £ — n E(n £ , n) 



dx 



n £ u £ ■ d t V x V £ ,8 + -n £ d t \\7 x ^ £ . s \ z - -A# M 
n \ 2 e 



»a,t. P(n £ )-P(n) 



dx 
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/ -n £ V x ^ £> s • (U £l5 -u e ) dx- \ f n £ \7 x $ £ ■ U E ,s dx. 
Jn t e Jn 



Moreover, we compute 

n E u E ■ d t V x ^^ s dx = p'(n) I N £ d t s £ .s dx - p'{n) 



t=s 



t=0 



-n / N E d t Ax 1 [s E<s ] dx + n 



N s A- l [s s> s] dx 



1 



t=o ' Jn 



N £ s £i $ dx 
n £ u £ ■ V x * e .5 da; 



and 



/ ^n e d t \W x ^ e , s \ 2 dx = f -(n e -n)d t \V x ^ e ,s\ 2 dx + f ^-nd t \\7 X V £ . S \ 2 dx, 
Jn 1 Jn 1 Jn 1 

where, by virtue of (fOjl . (f7T4|) . 
1 



-nd t \V x V £ ,s\ 2 dx 
n 1 

Next, we have 



P'{n)^d t \s £i s\ 2 dx+ — 
2 e 



V X * E , 4 • V^A" 1 ^-] dx-- \V x ^ £ j\ 2 dx 



t Jn 



- / A* 



E.S- 



P(n £ ) - Pin) 



dx 



p'(n)A^ £ . s N £ dx + 



P(n £ ) - ^^-{n £ - n) - P(n) ) Af £ s dx, 
n 1 



where 



- / Mf e , s N e dx = - / d t s £ , s N £ dx. 



Finally, 



n £ V x § £ ■ U £ , 5 dx = I iV £ V 
n 

£ 



■ U £ , s dx + - I \7 X 



V x ^ £ ,s dx 



J^N £ V X (^j - V £tS dx + n J^N £ d t A^s^s] dx. 



Summing up the previous considerations we may infer that 
jf f ~n £ |u £ - v V^e,^ 2 + ^E(n £ ,n) + ~ 



(s, •) dx 



(7.7) 



77 2 



+ / (p , (n)-\ S£ j\ 2 -p , (n)N £S£ , s + nN £ A] v 1 [s £iS } + — (-A„)~ 1//2 [s £ ,s] )(s,-)dx 



2 Jo 



V x (u £ - U £ , 5 ) + V^(u £ - XJ £iS ) - -div x (u £ - U £j(5 )I 



< 



2 1 

n - H[u ] - V x [*q,s]« + —E(n 0y£ , n 



dx dt 
1 



1 



o Jn 



n £ |u £ — U e j| dx dt 



dx+ 
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+ / (p'(n)-|[^o,e]5| 2 -p'(n)^o, £ [Afo,e]5 + niV 0!e A w 1 [iV ^] 5 + — (- A^ 1 / 2 [N , e ] s 



dx 



+ / Qe,S dt, 



where 



Qe,S 



[ n £ (d t v + u £ V x v + -v) • (U M - u e ) da; + / (n e d t V x ^ e ,s ■ v + n e u e ■ V x * e ,a ■ (U e ,$ - u e )) da; 
Jn v r y in 

2 e - n)U £)l5 • V x * e . 5 dx + - [ n e (u E - U El<5 ) • V^.a da; + ^ / (n £ - n)d t \V x t> e . s \ 2 dx 

T Jn * Jn 



+ / § E {V x U E ,s) : V X (U M - Ue) dx - / A* £ ,i 



j/(n) 



3 P K) 



P(n s ) - P(n) n s - n E(n e , n) 



4>, 



dx 



(n E - n) - P(n) A* e , s dx - / iV e V x — ■ U e , 4 dx. 



7.1 Dispersive estimates of the oscillatory component 

Our goal is to show that solutions s £i s, ^e,s of the homogeneous "acoustic" equation (|7.3|) . (|7.4p decay 
to zero in the L°° norm as e — > for any positive time t. To this end, we start with the total energy 
balance 



d 
dt 



[ (\V x * e , s \ 2 +p'(n)\s SiS \ 2 +n (-Aat)- 1 ^^] 2 \ dx + - [ \\7 x V e . s \ 2 dx = (7.8) 
Jn \ / T Jn 



yielding, in particular, existence and uniqueness of (weak) solutions to problem (17.31) . (|7.4[) provided the 
initial data are smooth and decay sufficiently fast for |x| — > oo. 

Taking advantage of the special geometry of waveguides, we consider the functions w k (z), z € B - the 
eigenfunctions of the Neumann Laplacian — Ajv,b in the (bounded) domain B C R 3 ~ L : 

-A N . B w k = A fe w fe in B, V z w k ■ n z \ dB = 0, A = < Ai < . . . < X k < . . . , k = 0, 1, . . . . 

The smoothing operators [g]s, g — g{x), x — [y, z] are defined as 

[9}s(y,z)= K s{v)* yP&{y) A k\g]{y))w k {z), (7.9) 



0<fe<l/<5 



where 



Mg}(y) 



\b\ 



ipg € C^°(R L ) is a cut-off function, 



o < i> s < 1, My) 



g{y,z)w k {z) dz, 

1 for \y\ < 1/5, 
for \y\ > 2/5, 
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and k$ is a family of standard regularizing kernels in the y— variable. 
A short inspection of (|T. 3[) . (17.41) yields 



s e ,s(t, x) = exp Se,S f~> x ) > ( 7 - 10 ) 

where s £ j is the unique solution of the Klein-Gordon equation 

df t s E: s -p'(n)As et s + ** e > s = °' ^ xSe - s ' n \ an = ®> 

emanating from the initial data 

5 e ,j(0, •) = = [M>, e ]j, ^S e , 5 (0, •) = -A[*o,J«. (7.H) 
n 

Consequently, thanks to the specific choice of the smoothing operators (17.91) . solutions s £i a take the 
form 

S e ,s(t,x)= ^ S k,e,s(t,y)w k (z), 
0<k<l/S 

where S k (t,-) solve the Klein- Gordon equation 

dl t S k , e , s - p'(n)A y S k , e , 5 + (n - ~ + p'(n)X k ^j S k , s>s = (7.12) 

for y belonging to the "flat" space R L , and with the initial data uniquely determined through (|7.11[) . 
Thus, employing the standard L 1 — L°° estimates for the Klein-Gordon equation 17.121 (see for instance 
Lesky and Racke [20l Lemma 2.4]), we have 

no u \\\ / <k) ( na , n , ao ,„ ,„ 





i + 3" 






2 



Going back to (|7.10[) we may infer that 

IK«(V)IU»(n) <u(to,e,S), te [to, T|, to > 0, (7.13) 

and, using (|7.3[) . 

||A* e ,*(t,-)IU~(n) <u(to,e,S), t£[t ,T], t > 0, (7.14) 

where w(to, £, <5) — » if £ — > for any fixed to > 0, 6 > 0. 
Finally, we claim the standard energy bounds 

||V x *e,5||z,~(o,T ; w™. 2 (n;ii3)) < c(m,5), (7-15) 

l|s£,<5lU~(o,T;i4»". 2 (n;_R3)) < c(m, S), m = 0, 1, . . . , (7.16) 
where the constants are independent of e for any fixed 5 > 0. 
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7.2 Asymptotic limit e — > 

Our next goal is to let e — > in (|7.7p . and, in particular, in the remainder Q e< s- 
1. We have 

y n e (d t v + u e V K v + ~v) • (U £>5 - Ug) da; 
= J n e (d t -v + U e , 5 V s v + ~vj ■ (U e ,a - u s ) da; + y n e (u e - U E)(5 ) • V K v • (U ei<5 - u g ) dx, 



where 



n e ( <9 t v + U £i(5 V x v + -v ) • (U e ,a - u e ) da; 



n E V x Il ■ (u e - U E ,s) dx + / n £ V x ^e,5 ■ V x v • (\J e ,s - u s ) dx, 



and 



/ n e V x II-(u e -U Ei5 ) da; = / n E V a II-u e dx-e / iV e V a; II-(v+V a; * e ,5) dx+n / VjII-V^^ dx. 
Jsi jo is] Jn 

Since 

where, by virtue of estimates (|4.1[) . (|4.3j) . 

[V^] 0SS% Ai7u £ -> SU weakly-(*) in L°°(0, T; L 2 (Q; R 3 )), div x U = 0, 

while 



[V^]r es V^u £ -> in L°°(0,r;L 5 / 4 (r!)), 



we get 



ess sup 

t6(0,Ti o 



leV^II • u e dx 



for e -> 0. 



Similarly, we use (|7.15l) to observe that 

jf H- J W X U ■ V x V e ,S dxj ^ in L 2 (0, T loc ) as e -> 



for any fixed (5 > 0. 
Thus we conclude that 



where 



j n e (d t v + u e V x v + -^v) ■ (U e>5 - Ug) dx < n £ |u e - U M | 2 dx + h\ 

hl s -> in i 2 (0, Tioc) as e -> 0. 



<■>■ 



(7.17) 
(7.18) 
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2. Taking advantage of the fact that div^v = we can write 

/ n e d t V x ^f e ,s ■ v da; = e / N e d t V x ^ e ,s • v Ax, 
Jn Jn 



where edt^x^e.s can be expressed by means of equation (|7.4[) . Using (IT. 15[) . f|T. 16[) we conclude 
that 



n e d t V x ^e.8 ■ v dx 



h 2 



with 



h 2 e 8 -> in £ 2 (0, T Xoc ) as e -> 0. 
3. Using ([7TT3]) . (I7TT5]) , we show that 

/ n e u e • V 2 ,*^ • (u E - U El<5 ) da; + / n £ (u £ - U e , 5 ) ■ 
Jo Jo 

+ / (n.-SjU^-VA^di + I {n e -n)d t \V x ^ EtS \ 2 Ax 
Jn Jn 

in L 2 (0, Tioc) as e 0. 



da; 



< ^ 



with 



4. Now, 



/ S £ (V X U E , 5 ) ; V»(U e ,a - u e ) da; 

Jf2 



< 



V x (u e - U M ) + V^(u e - U M ) - -di V;E (u e - U e , f )I 



da; 



+c/i e / |V X U £ . 5 | 2 dx. 



5. Next, in accordance with flU} and (1TT1)) , (|7T31) , 

f P(n e ) - P(n)n e -n E(n e ,n) 



a V 



d.r 



s- 2 , 



P( ne ) _ ?-^-{n £ - n) - P(n) ] dx 



(7.19) 
(7.20) 

(7.21) 



(7.22) 
(7.23) 



(7.24) 



< hi s -> m i 2 (0,T loc ) as e 0. 



6. Finally, using ([539"]) . ([Oj) . and ([jO]) . we infer that 



5 v a; , 



■ Ve,S Ax 



= h b eS ^0in 1/(0, T) as e -> 0. 



(7.25) 
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Using estimates (|7.17l - 17.25)) in (|7.70 we conclude that 

1 



n e |u e -v- V^^^l 2 (a,-) dx 
Jo \^ / 

- / ( -^■£ , (n e ,7i) -p'{n)N e s s ,s + p' (n)^\s e . s \ 2 ) (s, •) dx 



a V 2 



+ niV e A- 1 [ Se , a ] + y (-A n )-V2[ S£i5 ] ( Si .) dx 



< Qri , E |u , s - H[u ] - V x [* , e ]<s| 2 ) dx 
+ / ( ^E(no,e,n)-p'(n)N ,s[No,e}s+p(n)?- [N . E ] S 



dx 



n 2 



+c 



where 

Now, we claim that 



+ n7Vo, e A^[7Vo, £ ] d - + y (-A„) -1 / 2 [vYo, e ]a 
J ^ Qn e |u e -v- V,* e , 5 | 2 ^ dxdt + J^ h 6 £>s dt 
h 6 e S -> in L 2 (0, T loc ) as e -> 0. 



d.r 



ess sup 
te(o,T) 



E(n e ,n) , /y(n) 
2n 



as e -> 0. 



L 4 /3(n) 



Indeed we have 



E(n s ,n) p'(n) 



(H{n £ ) - H'(n)(n e - n) - ff(n) N \ 1/2 f 1 „„ (n e - n) 2 N 1/2 ' 

H (n) 



+ 



V e 2 y \2 e 2 

ff(n £ ) - H'(n)(n £ -n) - H(n)\ 1/2 (\ . An, -n) ^ 1/2 ' 



where 



^ H(n E ) - H'(n)(n E -n) - H(n) ^ 1/2 



lg"(n) ( " g -" )2 NJ 2 



while 



/ \2 ' ' 

in ^(O.T;! 9 ^)), 1 < q < 5/3, 



r 






ess 



(7.26) 



(7.27) 
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<V~e 



for a certain £ g [^/2, 2n]. 



Consequently, relation (|7.26p can be written in the form 

J Q« s |u e - v - V,* e , 5 | 2 ^ (s, •) da; 

E( "-- 77 ' ■ / ^ S£ , 5 V(^-)^+/i (-A)-^[^-iV e ] 2 ( S ,.)d, 



< Qri , E |u , s - H[u ] - V x [* , e ]<s| 2 ) 
+ / [-a£?(no, e ,n)-p'(n)JV , e [JVo J Jff [Md jE ]« 



















-he 



dx 



$0,< 



Jf! 



+ n7V , £ A 7V 1 [iV , £ ] 5 + — (-A n )-^ 2 [N , e ] s 
in e |u e -v- V^.al 2 ) dxdi+ / h 7 s At 



dx 



with 



h 7 e S -> in i 2 (0,Ti oc ) as e -> 0. 
Applying Gronwall's lemma we therefore get 

'1 



-n £ |u £ - v - V x ^ £:5 | 2 (s, •) da; 



J?(n.,») .M«)n \ (-> . )daj+ / 1 ( _ A) _ 1/2 ( _ ^ 2 ^ . } dx 



2 - ' v n 2 

-Ct jj \ „Cs 



< + C / i? £ , 5 (i)e- ct di e° s , 



o 



where 



jf Qn , e |u , s - H[u ] - V x [* , e ]<5| 2 J dx 



n 



+ / [ ^E(n , e ,n) - P '(n)N , e [N .e}s + p'(n)- [N , e ]s 



+ nN , e A N 1 [No, e ] s + — (-A n )-^ 2 [N 0je ] s 



d.r 



dx + I hl s dt 



(7.28) 



Thus, letting e — > we obtain 

esssup sup I (^-n e \u E -v-\7 x ^ eS \ 2 ) (s,-) dx <x(S), 
£^0 sS (0,Ti oc ) Jn \ z J 

where the function \ is determined in terms of the initial data, and x(^) — > as £ — > 0. 



(7.29) 
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7.3 Asymptotic limit <5 — ^ 



Letting S — » in (|7.29[) we may infer that 



H[V^ e Ue] -> VWv in L 2 (0, Ti oc ; L 2 (if; i? 3 )), 



(7.30) 



H- L [Vn e u e ]->0in L 2 (0, Ti oc ; L 2 {K ; i? 3 )) 



(7.31) 



for any compact if C f2. 

Relations (fPUl) . (|7T3T|) complete the proof of Theorem l3~2l 
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